Energy dependence of the dijet imbalance in Pb-Pb collisions at 2.76 ATeV 



Department of Physics, P. O. 
Helsinki Institute of Physics, 



Thorsten Ren^ 
Box 35, FI-40014 University of Jyvaskyla, Finland and 
P.O. Box 64, FI-OOOI4 University of Helsinki, Finland 



O 
(N 

< 
(N 



Oh 
I 

Oh. 
CD ■ 
^ ■ 



> 

^' 

o 



X 



The appearance of monojets is among the most striking signature of jet quenching in the context of 
ultrarelativistic heavy-ion collisions at the Large Hadron Collider. Experimentally, the asymmetry 
of back-to-back jets is quantified in terms of the dijet imbalance distribution by the ATLAS and 
CMS collaborations. Recently, the CMS collaboration has also studied the trigger jet momentum 
(Pt) dependence of the imbalance in the range between 120 and 500 GeV which is found to decrease 
with jet Pt. In this work, results from the in-medium shower code YaJEM are compared with this 
data set. These results suggest that the main effects observed in the data are the kinematical 
collimation of jets and the increase in the probability to produce more collimated quark jets with 
jet Pt, whereas there is no indication for any non-trivial energy dependence of the shower-medium 
interaction mechanism itself. The data furthermore can rule out models in which the jet shape is 
collimated due to the medium modification. 
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I. INTRODUCTION 



Hard partonic reactions in perturbative Quantum 
Chromodynamics (pQCD) lead (to leading order in the 
strong coupling Us) to highly virtual back-to-back parton 
configurations which subsequently evolve as parton shovi^- 
ers, hadronize and are then experimentally observed as 
back-to-back jet events. If such a hard process takes place 
in a medium as created in ultrarelativistic heavy-ion col- 
lisions, the evolving shower interacts with the medium 
and thus the properties of the jet are altered. 

Dijet imbalance observations as have been made by 
the ATLAS [l| and CMS [3 collaborations at the Large 
Hadron CoUider (LHC) for 2.76 ATeV Pb-Pb colhsions 
are an attempt to quantify the medium modification of 
the shower and several theoretical computations have 
been confronted with this data so far (see e.g. 
Since energy-momentum conservation guarantees that 
the shower-initiating partons are balanced in momen- 
tum (up to intrinsic fcy generated in the initial state par- 
tonic wave function and hard gluon radiation leading to 
three-jet topologies), the momentum imbalance is chiefly 
a statement about how much of the original parton mo- 
mentum is reconstructed given a particular jet definition. 
This in turn means that event-by-event fluctuations of 
the observed energy-momentum flow in a jet determine 
the imbalance distribution in a crucial way. 

Such fluctuations can have different sources: 1) fluctu- 
ation in the detector response 2) dynamical fluctuations 
of the parton shower pattern in vacuum 3) fluctuations of 
the soft background clustered into the jet in the presence 
of a medium and 4) modiflcations of the shower evolu- 
tion by the medium. 1) and 2) lead to sizeable dijet 
imbalances even in p-p collisions, and 3) increases this 
imbalance to some degree even in the absence of flnal 
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state shower-medium interactions. Any computation of 
the final state interaction requires hence a reasonable ac- 
count of the fluctuations in the baseline 1) - 3) , as back- 
ground fluctuations are potentially capable of masking 
the desired signal Q • The presence of background fluc- 
tuation is part of the reason for the lack of tomographic 
sensitivity of the set of flrst imbalance measurements as 
discussed in 

The CMS collaboration has subsequently also mea- 
sured the Pt dependence of of the imbalance distribution. 
This data has the potential to overcome the lack of tomo- 
graphic sensitivity mentioned above. The reason is that 
the different sources of fluctuations are expected to have a 
different Pt dependence — while detector response para- 
metrically scales AP^*^* ^ \/Pt, the background medium 
fluctuations occur at a constant scale AP!^^ . The dy- 
namical fluctuations due to the vacuum radiation have 
no simple form and depend on the jet deflnition. The 
main physics effect is the increased transverse collimation 
of jets at higher Pt both due to the kinematic forward 
boost and due to the increased probability of producing 
(narrower) quark jets, thus for given jet cone radius the 
fraction of jet momentum reconstructed inside the cone 
is expected to increase with Pt. Finally, the Pt depen- 
dence of the medium modiflcation depends on the precise 
scenario assumed for shower-medium interaction and is 
thus in principle tested by the data. 

The aim of this work is to illustrate these points using 
the example of the in-medium shower evolution Monte- 
Carlo (MC) code YaJEM [13, El- This code can be run 
simulating different physics scenarios, here we compare a 
multi-observable constrained scenario dominated by ra- 
diative energy loss with a small component of elastic en- 
ergy transfer into the medium (YaJEM-DE) j2j where jet 
shape is widened due to the interaction with the medium 
with a pure elastic drag scenario (YaJEM-E) [ll| in which 
the medium modification coUimates the jet shape which 
leads to a different prediction for medium-induced out of 
cone fluctuations. After a short overview over the model 
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framework, results for the energy dependence of the di- 
jet imbalance are presented and compared with the CMS 
data. This is followed by a discussion of the findings. 



II. THE MODEL 

The dijet imbalance is computed event by event us- 
ing a MC strategy. The framework is described in de- 
tail in [tI where it is applied to the ATLAS dijet im- 
balance measurement, here we summarize the essential 
points and highlight the physics relevant to understand 
the results. The basic outline of the modelling procedure 
is as follows: We first generate a hard back-to-back par- 
ton configuration based on a leading order (LO) pQCD 
code il2j, sampling both parton type and momentum. 
To effectively account for initial state effects and next- 
to-leading order effects, an intrinsic kx imbalance with 
a Gaussian distribution is introduced on the parton level 
(i.e. partons are never exactly back-to-back). The par- 
tonic configuration is, distributed with the transverse bi- 
nary collision profile, embdedded into a hydrodynamical 
medium which is constrained by reproducing LHC bulk 
multiplicity and spectra fis'l (note however that the re- 
sults of 7] strongly suggest that the dijet imbalance is 
insensitive to even gross features of the medium model 
beyond the mean density). 

Based on the parton path through the evolving 
medium, we next utilize the code YaJEM [13, Ull to com- 
pute the in-medium parton shower. The parton-medium 
interaction in YaJEM is cast into the form of transport 
coefficients at the parton position C which are q{C) (mean 
virtuality transfer to a parton per unit pathlength) and 
e(^) (mean energy transfer from parton to the medium 
per unit pathlength) . These are assumed to relate to the 
hydrodynamical energy density e as 

q[e]{C) = Kr[Ke] ■ 2 • e(C)'/"(coshp(C) - sinhp(C) cos V). 

(1) 

where p is the transverse flow rapidity of matter at po- 
sition C and ip is the angle between flow and parton direc- 
tion. K and Ke are adjustible parameters determining 
the overall strength of the interaction between partons 
and the medium. Note that YaJEM utilizes these trans- 
port coefficients to explicitly change parton kinenatics 
in the evolving shower in contrast to other in-medium 
shower models [H, [l^ which change the probability dis- 
tributions for soft gluon emission rather than kinematics. 

The ratio of Kn and Ke determines the relative con- 
tributions of medium-induced radiation and elastic en- 
ergy loss into the medium. Here, we explore two sce- 
narios: The default scenario, YaJEM-DE utilizes Ke — 
0.1K,Kji — 0.9K, i.e. about a 10% contribution of elas- 
tic energy loss. This scenario is the result of a combined 
analysis of multiple RHIC and LHC single hadron and 
hard dihadron correlation observables [l^. It leads 
to a substantial broadening of jet shapes below ~ 4 



GeV (i.e. at a scale comparable with the medium tem- 
perature) and leaves the jet structure almost unmodi- 
fied at higher momenta (cf. the scenario RAD in 17]). 
We contrast this with a pure elastic drag model where 
Ke — K',Kf( = in which the medium does not in- 
duce any radiation but in which a substantial momentum 
flow exists from the shower into the medium. This sce- 
nario leads to a coUimation of jet shapes (cf. the scenario 
DRAG in [13). Note that this scenario can be ruled out 
by several other high Pt observables [§] , however we use 
it to illustrate it to show the sensitivity of the energy 
dependence of the dijet imbalance to the jet shape. 

K (and K') are fixed as follows: K is determined by the 
requirement that Raa of single inclusive charged hadrons 
in central 2.76 ATeV PbPb collisions agrees with the data 
around 60 GeV (about the kinematical region where the 
lowest energy jet measurements are done). This leads 
to a R^f^ = 0.48 ± 0.02 in fine with ATLAS measure- 
ments of Rc p ^] and is consistent with K determinations 
from RHIC data within the uncertainty of extrapolating 
a hydrodynamical model from RHIC to LHC [13"] . Since 
YaJEM-E is not consistent with other observables, K' 
is selected to give the same trigger rate (i.e. the same 
R^*^) as YaJEM-DE. 

In order to compare with the measurement, the par- 
tonic output of the MC code in terms of event records 
must then be converted into jets. Hadronization is done 
using the Lund model [l9| . Since current jet measure- 
ments [3, [2] are not corrected for effects like detector 
response fluctuations or background fluctuations, a rig- 
orous comparison of theory with the data would require 
the following steps: 1) generation of a bulk event using for 
instance hydrodynamics with event by event fluctuating 
initial conditions 2) generation of a hard event inside the 
same bulk event (since there is an observable correlation 
between medium fluctuation hotspots and hard vertices 
the hard and soft physics can not be modelled inde- 
pendently in a rigorous treatment) and evolution of the 
perturbative shower inside the medium 3) computation 
of the back-reaction of the shower to the medium to ac- 
count for correlated soft physics inside the jet area (for 
instance elastic recoil of medium constituents) 4) simu- 
lating the detector response to turn particle event record 
information into calorimeter tower information 5) clus- 
tering and analysis of the calorimeter-level (or particle 
flow level) event as in the experimental procedure. 

So far, no theoretical model has realized this full, nu- 
merically very involved workflow (note also that detector 
response simulation codes are usually not available out- 
side the experimental collaborations). Nevertheless, it 
is important to construct a meaningful baseline calcula- 
tion before flnal state effects are applied which at least 
contains qualitatively the phenomena affecting the ex- 
perimental measurement and quantitatively reproduces 
strength and distribution of the fluctuations. A compar- 
ison with the published baseline (here PYTHIA + HY- 
DJET in [^) suggests that fluctuations are significantly 
more complex than Gaussian, thus this is a non-trivial 
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FIG. 1: Dijet momentum ratio distribution as measured by CMS and as computed both in YaJEM-DE and YaJEM-E for 
different trigger jet momentum ranges in central 2.76 ATeV Pb-Pb collisions. Shown in two panels is also the away side parton 
channel decomposition and the effect of kinematical coUimation (see text), all compared with CMS data (note that jet 
definitions are not identical). 



task. For this investigation, a baseline which reproduces 
the fluctuations in the published experimental baseline is 
constructed as follows: A hadron contributes to the jet 
momentum if: a) the hadron is at an angle of i? = 0.3 
or less with the jet axis b) the hadron species is one of 
TT+,TT-,TT°,K+,K-,p,p,'y c) hadron Pt > 1 GeV. The 
resulting jet momentum gets a contribution from ran- 
dom background fluctuations of a Gaussian distribution 
centered at zero with 14 GeV width (note that this is 
different from what has been used for ATLAS R = 0.4 
jets in 0] — since the jet definitions of ATLAS and CMS 
are not the same, there is no reason for the numbers to 
agree). Following the procedure applied by CMS, all jets 
with an energy below 30 GeV are discarded. 

While these conditions are certainly not identical to 
the result of a PYTHIA + HYDJET event propagated 
through a GEANT detector simulation and clustered by 
the anti-fcy algorithm [20j . the important thing to note is 
that they reproduce the pattern of fluctuations in the re- 
construncted jet Pt baseline well. Thus, while no strictly 
rigorous comparison with data is possible, one has a rea- 
sonable baseline to study trends and scaling patterns of 
the medium modifications. 

To sumarize, fiuctuations of the reconstructed jet en- 



ergy within the model are created by the following mech- 
anisms: 1) intrinsic kx imbalance of the partonic ini- 
tial state 2) dynamical fluctuations by vacuum radiation 
leading to e.g. hadron production outside the jet cone 3) 
uncorrelated medium background energy fluctuations im- 
posed in a Gaussian approximation 4) medium induced 
additional broadening of the jet shape, leading to en- 
hanced out of cone and soft hadron production. Note 
that only the last point is characteristic for the shower- 
medium interaction model. 



III. RESULTS 

The main result of this work is shown in Fig. [T] where 
the momentum ratio distribution of associate jet momen- 
tum Pt2 and trigger jet momentum Pti is shown for 
various trigger momentum ranges between 120 and 260 
GeV. The general trend in both model calculations and 
the data is that jets become less imbalanced at higher Pt- 
While the constrained scenario YaJEM-DE is in reason- 
able agreement with the CMS data throughout the entire 
momentum range investigated here, the jet collimating 
scenario YaJEM-E tends to underpredict the asymmetry 
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seen in the data in the lower trigger momentum ranges. 
For higher trigger momenta the quahty of the description 
of the data becomes comparable in both scenarios. 

At first glance, this seems to contradict the findings of 
where no distinction between YaJEM-DE and YaJEM 
could be established. Note however the difference in cone 
radius R = 0.4 vs. R = 0.3 as used here — since the two 
scenarios differ in jet shape, it is expected that they do 
not respond in the same way to a change in cone radius. 

The non-trivial trigger Pti dependence of the two sce- 
narios supports the expectation that this set of data 
probes the momentum dependence of various fluctuation- 
generating mechanisms in a non-trivial way: Just based 
on the 260 - 300 GeV result, it would be difficult to dis- 
tinguish the two scenarios, however the full momentum 
dependence provides the information to overcome this 
problem. 

To further investigate the physics underlying the re- 
sults, for two trigger ranges results with the away side 
shower-initiating parton forced to a gluon or to a quark 
computed in YaJEM-DE are shown. These results are 
in line with the expectation that the broader transverse 
shape of gluon jets should lead to a stronger asymme- 
try by inducing more out of cone radiation. While at 
low trigger Pti the full result is somewhat closer to the 
gluon result, at high Pti this is no longer true and the 
full curve is about the average between quark and gluon 
result. This confirms the expectation that the rise in the 
relative fraction of quark jets with jet momentum con- 
tributes to the coUimation of observed jets and hence to 
a reduction in the observed imbalance. 

A different effect is the kinematical collimation of high 
Pt jets due to the forward boost of evolving shower. We 
can assess the role of this effect by forcing jets triggered 
at 220 GeV to have the same jet shape distribution as 
jets at 120 GeV. If this is done, as expected an increased 
imbalance is found which is of the same order as what is 
induced by the gluon to quark transition. We can con- 
clude from this that the decrease of the dijet imbalance 
with trigger momentum is caused equally by both effects. 



IV. CONCLUSIONS 

The results suggest that the observed change in the di- 
jet imbalance as a function of trigger jet momentum can 
be accounted for in a scenario that is well-constrained 
by other observables and can be explained in terms of 
essentially known physics, i.e. the increased collimation 
of jets due to kinematics and a transition to a less gluon- 
dominated regime. In other words, no particular spe- 
cial or exotic assumptions beyond kinematics need to be 
made for the Pt dependence of the shower modification 
by the medium itself. 

The failure of YaJEM-E to account for the data in 
part of the momentum range indicates that the jet Pt 
dependence of the imbalance is as expected sensitive to 
the jet shape and its distortion by the medium modifi- 
cation. Comparison with Q suggests that smaller cone 
sizes may offer observables more sensitive to the trans- 
verse structure generated by shower-medium interaction 
scenarios. 

In summary, the reasonable agreement between the 
trends of the imbalance evolution as seen in the data 
and in YaJEM-DE strengthen the case for the interpre- 
tation of the combined data of high Pt observables as 
being dominated by pQCD radiative energy loss with a 
small fraction of elastic energy transfer into the medium. 
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